function [x_next, phi_next] = SOIM(a, x, phi, n)
% SOIM 定义 SO-IM 映射的单步更新公式
% 说明: 此处 a 表示耦合参数 b
%
% x_next = 0.2 * (-1)^n + b * sin(pi * phi) * x
% phi_next = phi + x

    x_next = 0.2 * ((-1)^n) + a * sin(pi * phi) * x;
    phi_next = phi + x;
end
